Four-qubit device with mixed couplings 
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We present the first experimental results on a device with more than two superconducting qubits. 
The circuit consists of four three-junction flux qubits, with simultaneous ferro- and antiferromagnetic 
coupling implemented using shared Josephson junctions. Its response, which is dominated by the 
ground state, is characterized using low-frequency impedance measurement with a superconducting 
tank circuit coupled to the qubits. The results are found to be in excellent agreement with the 
quantum- mechanical predictions. 



PACS numbers: 85.25.Dq, 85.25.Cp, 03.67.Lx 

The implementation of one- and two-qubit gates in su- 
perconducting qubit prototypes [1] has confirmed their 
utility for quantum computation. Such qubits are read- 
ily fabricated, and highly scalable in principle. A promi- 
nent subcategory consists of superconducting loops inter- 
rupted by three Josephson junctions, so-called 3JJ flux 
qubits [2]. If threaded by a near-degenerate magnetic 
bias flux <£> x w ($0 is the flux quantum), such de- 
vices have quantum states which are superpositions of 
clockwise (| [)) and counterclockwise (| j)) circulating su- 
percurrent J p . The small area of 3JJ qubits reduces their 
coupling to environmental magnetic noise, making com- 
paratively long coherence times possible [3] , but also lim- 
its the strength of their inductive coupling [2, 4]. This 
can be overcome using direct galvanic coupling through 
a shared Josephson junction [5]. When reporting its ex- 
perimental realization [6] , we mentioned that direct cou- 
pling can have additional advantages. First, the coupling 
strength can be varied independently of the sample ge- 
ometry by changing the shared junction's critical cur- 
rent. Second, while direct coupling normally has the an- 
tiferromagnetic (AF) sign just as in the inductive case, a 
"twisted" design (joining two qubit loops in a "00" shape 
with crossing leads in the center, see Fig. lb) features fer- 
romagnetic (FM) coupling [7]. 

In this Letter, we pursue this by studying a four-qubit 
circuit in which the two types of coupling co-exist. This 
is very promising from the perspectives of realizing non- 
trivial Ising-spin systems [8] and scalable adiabatic quan- 
tum computing (AQC) [9, 10]. However, while our cir- 
cuit behaves in full agreement with quantum mechanics 
and features excellent multi-qubit bias control, the mea- 
surements presented are essentially equilibrium. The sys- 



tem's effective low-energy (pseudospin) Hamiltonian is 
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where is the bias on qubit i (implemented through a 
flux bias <£ x — |$o) 5 \ is its tunnelling amplitude, Jij is 
the coupling energy between qubits i and j, and a z , <j x 
are Pauli matrices in the span of \l) (a z = — 1) and ||) 
(a z = 1). Equation (1) adequately describes a system of 
flux qubits [2, 4, 11], although possible limitations [12] 
remain to be investigated. 

Our approach follows earlier work using the impedance 
measurement technique (IMT) [13]. This method mea- 
sures impedance shifts in a high-quality superconducting 
LC (tank) circuit inductively coupled to the system. For 
one qubit [14], the average loop current / = (|6| 2 — |a| 2 )/ p 
in an eigenstate a\l) + 6||) changes direction at the anti- 
crossing point <1> X = (where \a\ — |6|, i.e., the pseu- 
dospin ground state flips from ||) to ||) and vice versa for 
the excited state). This causes a peak in the qubit suscep- 
tibility x = dl/d§*; since dl/d$>* ex d 2 E/d(§*) 2 , where 
E is the average energy and <i$x is adiabatic (note 15 in 
Ref. [14]), this also follows from the large energy-band 
curvature at the anticrossing. Note that the very pres- 
ence of such peaks is already a quantum effect, since 
the classical x($ x ) curves feature hysteretic loops [10] 
in the absence of flux tunneling. Due to the qubit- 
tank coupling, x affects the tank's eigenfrequency cut and 
hence, under resonant driving, its current-voltage phase 
angle 6 ex x- Thus, the width and height of such an 
IMT-peak in 6 provide information about the anticross- 
ing, i.e., about A and 7 P . 
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FIG. 1: (a) Electron micrograph of the sample. The central 
junctions A1-A3 couple the Al qubits qi-q^. The surrounding 
Nb coil is part of the LC tank circuit, used for both mea- 
surement and global flux biasing. The Nb lines it>i-b4 allow 
asymmetric bias tuning, (b) Schematic diagram. The driving 
current idr(t) is much smaller than the oscillations in It and 
only needs to compensate for the latter 's losses. 



For multiple qubits, a more general analysis yields [15] 
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where Lt and Qt are the tank's inductance and qual- 
ity factor respectively, = e~ E »l T j S ^ JV e~ E »l T is the 
Boltzmann factor of the eigenstate of H in (1), and 
with the matrix element 
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Xi are tank-qubit coupling coefficients in H[ nt = 
I r f J2i ^i a z l \ and J r f is the oscillating component of the 
tank current It = IbT + Irf(t). The term i = j pro- 
duces the IMT-peak corresponding to qubit i. Cross 
terms i ^ j, on the other hand, appear only if or 
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FIG. 2: (a) -tan8(/ b T) at I hl = 700 jiA, I h2 = 10 /iA, and 
ib4 = 250 fiA (ib3 = throughout). The four qubit peaks 
do not overlap, (b) The same graph (vertically shifted for 
clarity) but at ibi = 268 /xA, where the peaks corresponding 
to qi and (73 overlap. The height increase compared to the 
sum of the two non- overlapping peaks is a signature of FM 
coupling. Dashed lines are theoretical fits. The fit for q 2 is 
poor because the relatively large tank amplitude washes out 
the experimental peak [13]; a better fit has been obtained for 
smaller amplitudes (not shown). 



\v) has nonzero entanglement between qubits i and j due 
to their interaction. For two weakly coupled qubits [4], 
the effect of interaction can be understood as a correc- 
tion. The cross terms in R^ v then change the shape of 
overlapping peaks at co-degeneracies [i.e., when the in 
(1) are chosen such that both qubits flip simultaneously], 
which turn out to be smaller (larger) than the sum of the 
two individual peaks for AF (FM) coupling [16]. 

Here, the analysis needs to be extended both beyond 
two qubits and to stronger coupling (particularly J > A), 
which qualitatively changes the locus of the IMT-peaks, 
as can already be understood classically: neglecting the 
Ai in (1), the degeneracy points of the classical multi- 
qubit flux states are seen to depend on the Jij [6]. One 
then proceeds by fitting the coefficients of H in (1), see 
below. 

Figure la shows the sample with four Al qubit loops 
inside a Nb pancake coil. The areas of the Al loops are 
10 x 7.5 jim 2 . Two junctions in each qubit are typically 
650 x 150 nm 2 , while the third is -25% smaller [2]. The 
Nb coil has Lt = 105 nH, and together with an external 
capacitance Ct = 470 pF forms a parallel tank circuit 
with uj t /2tt = 22.703 MHz and Q T = co T R T C T = 400, 
where Rt is the effective resistance. The Al loops 
were fabricated by e-beam lithography and conventional 
shadow evaporation, and the Nb coil by e-beam lithog- 
raphy and CF4 reactive-ion etching. 

The qubits are coupled to each other both magnet- 
ically and via a shared junction, consisting of parallel 
junctions A\—A^ (Fig. 1). For a coupling junction with 
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FIG. 3: (a) — tan 0(it>T, it>i) at I h2 = 10 /iA and I h4 = 
250 \i A. The repulsion between the traces (^1,^2), (^2,^3), 
and (#1,(74) shows AF coupling between the qubits in those 
pairs, while the merging of traces (^1,^3) demonstrates FM 
coupling. The line widths are proportional to the tunnelling 
amplitudes of the individual qubits. Since the qubits' sensi- 
tivity to ibi decreases in the order q±, (72, #3, q<± (cf. Fig. la), 
the various slopes allow one to identify the trace belonging to 
each qubit. (b) Theoretical graph for the same parameters. 



large critical current 7 C , and assuming identical qubits, 
J = hl*/2el c [6]. For our system, I p ~ 250 nA. Due to 
the relative twist [6] between the qubit loops (Fig. lb), 
the couplings qi~q3 and q2~q4 are FM (J^- < 0), while 
the others are AF (J^- > 0). The biases ^ can be varied 
by adjusting the currents ibi-b4 through the Nb wires, 
as well as the dc current ibT through the tank coil. The 
qubit-coil mutual inductances follow from the flux peri- 
odicity as M qliT ~ M^t ~ 38.2 pH, M q2 , T = 38.6 pH, 
and M g4 ,T = 38.9 pH. 

The sample was thermally anchored to the mixing 
chamber of a dilution refrigerator at a temperature 
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FIG. 4: Same as in Fig. 3, but as a function of ibT and ib2, 
at ibi = —400 \i A and ib4 = — 50 fiA. The skipped crossing 
at the center corresponds to four-qubit co- degeneracy. Due 
to their FM coupling, the qi- and (73-traces merge over an 
extended region. The qubits' sensitivity to ib2 decreases in 
the order #2, <?4, <?i,3- 

T m i x ^10 mK. In general the effective T is affected by 
noise on external leads and so is larger than T m i x . Of- 
ten, T m i x alone is given, since T is difficult to determine. 
One of the advantages of IMT is that it allows one to 
determine effective temperatures [6, 14]. From the best 
theoretical fit we estimate T « 70 mK, well within the 
expected range. 

Figure 2 shows peaks in — tan O as the overall flux is 
scanned by changing IbT- For most biases, four peaks can 
be distinguished (Fig. 2a), each corresponding to degen- 
eracy in one of the qubits while the others provide a semi- 
classical static background field. The narrow peak for q<i 
indicates a small A2 (~12 mK), and incidentally con- 
firms that the width of the other traces is not resolution- 
limited. When qi and qs are co-degenerate (Fig. 2b), the 
resulting peak exceeds the sum of their individual peaks, 
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indicating strong FM coupling as mentioned above (this 
is also true for q 2 and q±). 

Figure 3a plots — tan 0(ibT?^bi)- The four traces are 
IMT-peaks, each corresponding to a single-qubit anti- 
crossing. Since JbT biases all qubits almost equally, the 
trace with the greatest sensitivity to ibi can be ascribed 
to qi etc. Combining several such cross-sections of the 
total IMT response, one reconstructs the biasing coef- 
ficients in = M qiiT I hT + E^kb/bj- The cen- 
ters of the traces mark the boundaries of the classical 
stability diagram. Each region between the traces cor- 
responds to a different (minimum-potential) flux state 
of the qubits [17]. The shifts of the traces after each 
(skipped) crossing depend on the J^-. The peak shapes, 
on the other hand, carry information about the A^. As 
expected, the traces corresponding to two qubits with AF 
(FM) coupling repel (attract) each other, and in the FM 
case merge over a certain distance. Fitting the data to 
Eq. (1) yields Ai = 147, A 2 = 12, A 3 = 163, A 4 = 165, 
and J12 = J34 = 163, J14 = J23 155, J13 = J24 = -62 
(all in mK), and 
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where the missing entries correspond to Jb3 which was 
never used. The theory (Fig. 3b) accurately describes 
the system's behavior over the whole parameter space, 
especially at the co-degeneracies. 

Figure 4a shows — tanO(/bT, ^2) for a case when all 
four qubits become degenerate simultaneously in the fig- 
ure's center, where they produce two repelling peaks. 
From H in (1) one infers that, between those peaks, 
the two lowest eigenstates are close to superpositions of 
I HID and I HIT), for which all the coupling energies are 
negative; at co-degeneracy, the superpositions become 
(anti) symmetric [18]. From the viewpoint of AQC, it 
is encouraging that the ground state of H thus is glob- 
ally entangled [19]; for a quantification, see [20]. With 
all parameters already determined from Fig. 3 and a few 
other cross-sections like it, the theoretical comparison in 
Fig. 4b (and many other nontrivial cross-sections) can 
be regarded as having no free parameters, making the 
agreement remarkable. 

Complete agreement between experiment and theory 
was also obtained for a second sample, with Ai = A2 = 
60, A 3 = 130, A 4 = 110, J12 = J23 = 300, J14 = J34 = 
330, and J13 = J24 = — 90 (all in mK). However, the 
lowest effective temperature attained was only ~300 mK. 

In conclusion, ferro- and antiferromagnetic couplings 
between four 3JJ flux qubits have been realized simulta- 
neously with shared Josephson junctions, with a coupling 
strength significantly exceeding the inductive one. (Inci- 
dentally, the results also show that direct galvanic cou- 
pling can be used for more general circuits than the linear 



arrays considered in Ref. [21].) The data fully agree with 
a quantum-mechanical description to the experimental 
accuracy. Currently, quasi-equilibrium impedance mea- 
surement already provides valuable information comple- 
mentary to state readout, notably the determination of 
the Hamiltonian. With a faster and possibly qubit- 
selective detection method, the circuit used looks very 
promising for studying adiabatic bias manipulations. 
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